Pair correlations of a spin-imbalanced Fermi gas on two-leg ladders 
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We study the pair correlations of a spin-imbalanced two-leg ladder with attractive interactions, 
using the density matrix renormalization group method (DMRG). We identify regions in the phase 
diagram spanned by the chemical potential and the magnetic field that can harbor Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO)-like physics. Results for the pair structure factor, exhibiting multiple 
pairing wave-vectors, substantiate the presence of FFLO-like correlations. We further discuss phase 
separation scenarios induced by a harmonic trap, which differ from the case of isolated chains. 



The experimental realization of fermionic superfluids 
in ultracold atom gases under clean conditions and with 
a great control over interactions has paved the way to- 
ward a detailed understanding of the BEC-BCS crossover 
of spin-balanced, ultracold Fermi gases [1]. Now, the case 
of a polarized two-component Fermi gas, realized by un- 
equally populating the two lowest hyperfine states, has 
moved into the focus of current experimental work [2]. 

Intriguing properties such as phase separation in a 
trap and the transition from superfluidity to a normal 
state driven by the population imbalance have attracted a 
great deal of attention, but it is, in particular, the search 
for exotic superfluids such as the FFLO one [3, 4] that 
drives the current interest in imbalanccd Fermi gases. In 
an FFLO state, the order parameter is spatially inhomo- 
geneous with Cooper pairs with a finite center-of-mass 
momentum. In recent experiments on three dimensional 
(3D) ultracold gases, this state remains elusive, and the- 
oretical work indicates that in 3D, the phase space vol- 
ume of this phase in the interaction-polarization plane 
is small [5]. Yet, reducing the spatial dimension renders 
this pairing mechanism more effective as a larger portion 
of the Fermi surfaces of minority and majority spins can 
be matched [6]. In particular, in one dimension (ID), 
where a true condensation is prohibited, the existence of 
an FFLO-type state with quasi-long range order has been 
proven by means of analytical [7, 8] as well as numeri- 
cally exact approaches [9-12]. This also pertains to the 
experimentally relevant case of a harmonic trap [10-12]. 

It is then natural to ask whether such quasi-FFLO 
states seen in ID are stable against coupling chains to 
2D or 3D ensembles, in order to connect the aforemen- 
tioned theoretical results for ID to those available for 2D 
[6, 13-16]. For ID chains weakly coupled to a 3D array, 
the polarization-interaction phase diagram has recently 
been derived in Refs. [6, 7, 12, 17]. Here we present 
a rigorous and quasi-exact numerical analysis of pairing 
correlations in two coupled chains, using DMRG [18]. 

The physics of spin [19] and Hubbard ladders [20] has 
proven to be unique and interesting in itself, due to the 
emergence of exotic quantum phases - such as spin liq- 



uids - driven by strong correlations [19]. The experimen- 
tally realization of ladders in optical lattices as arrays of 
double wells could be an important step toward under- 
standing the experimental and theoretical challenges that 
wc may face before scaling to 2D systems. Double wells 
have already been studied in recent experiments [21] and 
Danshita et al. [22] have given a prescription of how the 
parameters of a double well potential need to be tuned 
to create a ladder system in an optical lattice. Methods 
to detect the FFLO state in experiments using noise cor- 
relations or interferometry are discussed in Refs. [12, 23]. 
We study the Hubbard model on a 2 x L ladder: 
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where c^^} ^ is a fermionic annihilation (creation) operator 
acting on a site on rung i and leg I {I = 1,2). The hop- 
ping matrix elements along rungs and legs are denoted by 
t_L and i|| , respectively. We study the case of strongly at- 
tractive onsite interactions U = — 8i|| (see Ref. [22] for a 
discussion of how U , ty, and i_L are related to experimen- 
tally controllable parameters in an optical lattice). Fur- 



ther, ni,i_, 



yielding the number of fermions 



of each species as = X]/ i{'ni,i,a), with N = + 
and the pseudo-spin index cr =t,i. As customary, h is 
the magnetic field, fj, the chemical potential, n = N/{2L) 
the filling factor, and p = (iV| — Ni)/N measures the 
polarization (we use > N^). A harmonic trapping 
potential is introduced by the last term in Eq. (1). 

Our analysis is mostly concerned with the case of t± = 
t|| and thus not restricted to weakly coupled chains. We 
first study the emergence of FFLO-like correlations as a 
fimction of the ratio t±/t\\ and filling n for Vtrap = 0. 
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FIG. 1: (Color online) Bulk phase diagrams {U = —Sty, t± — 
in the magnetic field h vs. chemical potential fi plane 
{L = 30 data): (a): filling n; (b) polarization p (see the 
legend for the color coding; dotted lines: guide to the eye). 
The phases I-V are defined in the text. Unlike the U > case 
[9] , no plateau at p oc hole-concentration is found. 



The bulk phase diagram is displayed in Figs. 1 (a) and 
(b), showing contour plots of n and p in the h vs. /i 
plane, respectively. We infer the presence of five phases 
(for the single chain case, see Rcf. [25]): I: the empty 
phase (n = 0, p = 0); II: a spin gapped phase {p = 0, n < 
1); III: fully polarized, but less than half filled {p = n, 
n < 1); IV: partially polarized (n < 1); V: fully polarized 
(n = 1). Phase IV is the candidate for FFLO-like physics 
and we will show that a rich structure in the real-space 
modulation of the pair correlation function emerges. 

As a guidance for the interpretation of our DMRG re- 
sults, we next consider the noninteracting limit (see, e.g., 
Rcf. [9]). The dispersion relation at [/ = with pe- 
riodic boundary conditions along the legs of the ladder 
is: e'^(fc||) = — 2t|| cos(fc|| ) — t ± cos{k'^) , i.e., there are 

0,7r. In 



four bands which we label with (fc J ,a); k° 
the imbalanced case, the Fermi surfaces arc mismatched, 
and we compute the pairing momenta Q = iQ\\,Q±) 
of possible FFLO-like states from the difference of the 
Fermi- wave vectors for majority and minority spins, i.e., 
Q — ^ = fc^ii — and Q± = kp^ ^f±' 

Depending on p, n, and t±/t^\ , the band structure can al- 
low for more than one pairing wave vector Q||, yielding a 
quasi-condensate with multiple contributing modes. This 
may also be thought of as coexisting quasi-condensates 
with Q± = and Qj_ = tt. Examples for the expected 
polarization dependence of the pairing momenta Qn are 
displayed with lines in Fig. 2 for generic fillings of n < 1 
and 1 1 = 1 and 2. To label these branches, we use a 

* ) 

^ , which specifies as well as the k'J_ of 

the majority and the minority spins. For instance, tt^'^'*^' 

represents Q± = tt pairs, formed from majority spins 

t .... I 

with kj_ = TT and minority spins with kj_ = 0. 

The disappearance of certain branches is due to a band 
getting completely filled or depleted, causing kinks in the 
other branches at the same p. In the case of the ladder, 
there are two constraints on the four Fermi momenta. 



symbol Q\ ^ 



71/2 



n/2 



_(a) t^=t||, n=0.4 

□ Q^=0 

(o.n) n „ 
It y^-it 

(n.O) 

O*"-"' ^.rt-^^ 

-IVb Ji^J^°IVc ^ 


(b) t^=t||, n=0.7 


-(c) t^=t||, n=0.9 

IVd / 

1 .""T^w*" 1 ,1,1, 


(d) t^=2t||, n=0.7 

rfjf" IVd ' 



)0' 



0.2 0.4 0.6 0.8 D 0.2 0.4 0.6 0.8 1 
P P 



FIG. 2: (Color online) Pairing wave-vectors extracted 
from (i) Q|| — fej^n — k^^^ at (7 = (lines) and (ii) the po- 
sition of the maximum in the pair structure factor (DMRG 
for L = 50 at f/ = — Sty; squares: Qi_ = 0; circles: Q± — tt). 
(a)-(c) t|| = t± with n = 0.4,0.7, and 0.9, respectively, (d) 
t± — 2t||, 71 = 0.7. The lines [see the legend in (a)] represent 
the possible branches (see the text for details). 



Luttinger's theorem for both spin flavors: X^fc^ \\ ^ 
2tt Nc/{2L) [20, 24]. Thus, unlike the single-chain case 
where n and p determine the two Fermi momenta inde- 
pendently of U , in the case of a ladder, interactions may 
modify the derived from the ?7 = case. 

Let us now come to the discussion of the (s-wave) pair- 
ing correlations at \U\ > 0. To that end, we compute 
the pair structure factor Fourier transform of 



pair 
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2L ^ 
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"hh,r 
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h,l2 



This quantity is displayed in Figs. 3(a) and (c) for k± = 
0, TT, respectively (n = 0.7 and t± = i||). For k± = Q [see 
Fig. 3(a)], there is a strong coherence peak at = at 
p = 0, which, upon polarizing, shifts to finite momenta 
(5|| > 0, as expected for FFLO-like pairing. In the k± ~ tt 
channel, a peak with Q\\ > emerges only for p > 0.2, 
approaches Q\i = at p w 0.34, and then increases 
again in the limit of large polarizations. 

We can further elucidate this behavior by analyz- 
ing the single-particle momentum distribution function 
(MDF) n^ii , which is the Fourier transform of the one- 

J. We dis- 



particle density matrix p^^^^ ■■ = (cj 



z,/i,cr j.l: 

play the MDF, for the parameters of Figs. 3(a) and (c), 
in Figs. 3(b) and (d). Most notably, we see how the frac- 
tion of majority fermions increases for both k± = and 
k± ~ TT, with the k± = channel dominating. Secondly, 
as expected, the minority spins in the k± = tt channel are 
depleted first at around p > 2/7. Therefore, for instance, 
a 0('^''^) quasi-condensate can only be realized at small p. 

In principle, the case of Fig. 3 could allow for up to four 
combinations of {Qh, Qi_) at small p, as can be inferred 
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FIG. 3; (Color online) t± = t\\, U = -8t\\, and n = 0.7: 

(a) , (c): Pair structure factor for (a): k± = (c): k± = tt. 

(b) , (d): MDF nl^^^,,^ {(^ =T(i): solid(daslied) lines). In all 
plots, p = 0,1/7, ... , 6/7 (except for a —[ in (d); there, p — 
0, 1/7, 2/7, 3/7). Arrows indicate increasing p. The legend in 

(c) applies to (a),(c), the thick solid line in (b),(d) is for p = 0. 



from Fig. 2(b). Yet, at J7 < 0, wc find at maximum one 
peak in the pair structure factor in each channel. We 
compare the position of these peaks as extracted from 
the DMRG data to the Q\\ computed in the U ~ limit, 
for the parameters of Fig. 2. In the case of t_L = we 
choose the fillings such that several cuts through phase 
IV of Fig. 1 can be followed (as indicated by the labels 
IVa-IVd in Fig. 2). We subdivide phase IV according to 
how many bands are partially filled (i.e., < n < 1): IVa: 
four bands; IVb: (0, a =1, |) partially filled [(tt, a =t, i) 
empty]; IVc: three bands [(tt, |) empty]; IVd: two bands 
[(7r,T) and (0,i); (tt, i) empty; (0,T) full]. 

At ?? = 0.7 [Fig. 2(b)], the Q± = quasi-condensate at 
small p starts on the O'^'^''^^ branch and then continously 
moves to the O*-"'*^^ branch. At larger p, a second con- 
densate with Q± ~ TT builds up (circles), namely one of 
the n^'^'^^ type. Again, the positions Q|| of the quasi- 
condensates well agree with the U ~ predictions (lines 
in Fig. 2) at polarizations p > 0.3. In the case of a large 
fiUing n ~ 0.9 and t± = t\i [see Fig. 2(c)], the picture is 
very similar to the n ~ 0.7 case, with the difference that 
at large p, one enters into region IVd, where the (0, t) 
band is filled, and hence only the n^^'^^ branch survives. 

A simpler structure for the quasi-condensate is ex- 
pected for the parameters of both Figs. 2(a) and (d), 
where at maximum, two modes may exist. In these two 
cases all possible (Qy , Q±) indeed contribute to the quasi- 
condensate at finite U and, moreover, their momenta Q|| 
mostly agree with the U = predictions. An interesting 
effect occurs at polarizations p > 0.35 in Fig. 2(d), where 
atU = 0, no phase IVc is present. The interactions renor- 
malize the dispersion such that the (0, t) band remains 
partially filled up to p < 0.7, allowing the O'-*''"^ branch to 
survive up to that value of polarization. Therefore, due 
to the interactions, a small window for the emergence 



of IVc opens up. For t± > 2t\\ and n = 0.7, only one 
quasi-condensate is realized (not shown in the figures). 
We emphasize that at any polarization and in all cases 
studied, at least one mode with Qn > is present. 

At small p, the decay of pair correlations is consistent 
with \pI%\ ~ |cos(Q|||r - s|)|/|r - s|« for p < 2/7, as 
expected for a single-Q quasi-condensate [7]. The same 
applies to all cases in which a single branch is present in 
Fig. 2. In all other cases, our data suggest that the en- 
velope of Iph'^^^I follows a power-law. Moreover, we find 
that pair correlations decay slower than those of compet- 
ing instabihties such as density-density correlations [12] 
at small p, while at present, we cannot make a definite 
statement about larger polarizations. Out of the possi- 
ble quasi-condensates, typically the ones with the small- 
est Q|| win and show up with a significant weight in the 
pair structure factor. The most ID-like behavior with a 
single {Q\\,Q± = 0) mode is encountered in phase IVb, 
which for t±/t\\ — 1,2 extends up to full polarization for 
n < 0.4. In IVb, Q|| = 2Tmp, as in ID [10]. 

We now turn to the effect of a harmonic trapping po- 
tential on the pairing correlations and the density profile 
by setting Vtrap = 0.002t|| with = 40 fermions and 
tj_ = As for phase separation, our results displayed 
in Fig. 4 indicate that at small polarizations {p < 0.2), 
there are three shells: an inner core with a vanishing 
polarization, a thin shell that is partially polarized with 
increased spin fiuctuations, and fully polarized wings [see 
Fig. 4(b), p = 0.1 there]. This is distinctly different from 
the case of uncoupled chains with the same U and filling 
[8, 10], where FFLO-like correlations develop in the core 
as soon as p > 0. By invoking the local density approxi- 
mation, we trace this back to the slope of the boundary 
between phases II and IV that has a different sign in the 
case of ladders (see Fig. 1) as compared to chains [25]; 
thus the system goes through IV in the center of the trap 
before entering III toward the edges. At U = — 8f||, the 
formation of fully polarized wings sets in as soon as p > 0. 

At larger polarizations, there are only two phases: 
in the core, the local polarization pi^i = (ri;,;,! — 
'^i,',i)/("i,',T increases and the unpolarized phase 

gives room to a partially polarized one, while the wings 
remain fully polarized [see Fig. 4(c) for the example of 
p = 0.6]. We relate the phases emerging in the trapped 
situation to those of the bulk system shown in Fig. 1: 
the unpolarized phase is a Qy = superfiuid - phase II 
- while the partially polarized one has, similar to phase 
IV, FFLO-like features. Our results for the pair struc- 
ture factor, presented in Fig. 4(a), confirm this picture: 
For p < 0.2, a coherence peak at Qn = dominates in 
the fcj^ = channel, which, as soon as the core assumes a 
finite polarization, develops into a > peak, charac- 
teristic of an FFLO state. The evolution of the coherence 
peaks for both k± ~ and tt resembles that of the un- 
trapped system [Figs. 2(a) and (b)]. The coherence peaks 
at both Q± ~ and tt survive up to saturation. At high 
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FIG. 4: (Color online) Harmonic trap Vtrap ~ 0.002t||; — 
t±, N = 40, and U — —Sty. (a) Pair structure factor for 
k± — and p = 0, 0.1, . . . , 0.9. (b), (c): Density profile 
{ni,i) {uij = ni,i,| + riij^i), local polarization and lo- 
cal {{slif) = - ni,,,i)2>/4 for (b) p = 0.1 and (c) 
p = 0.6. Vertical, dotted lines are a guide to the eye. 



densities (not shown here), a band- insulating Fock state 
with uncorrelated tightly bound pairs appears in the cen- 
ter of the trap, displacing the other phases toward the 
edges, similar to the ID [10] and the 2D cases [13]. 

Discussion - The physics of ladders already features 
some characteristics of higher dimensional systems. In 
particular, the positive slope of the boundary of phase II 
(see Fig. 1) is a remarkable effect of the dimensionality 
that changes the order of the concentric phase-separated 
shells in a trap, in agreement with mean-field calculations 
in 2D and 3D. To learn about the incipient 2D physics, 
it is edifying to look at a ladder with t± ~ 2i||. In this 
case, the IVb and IVd regions occupy most of phase IV, 
with only a small sliver of IVc at high densities due to the 
renormalization of the bandwidth [see Fig. 2(d)]. In the 
isotropic 2D system, region IVb would correspond to two 
partially filled bands, while IVd has a band of majority 
spins above half filling, touching the boundaries of the 
Brillouin zone. At a low density (phase IVb), the physics 
fits within a single-band picture [see, e.g., Fig. 2(a)], sim- 
ilar to a chain. At larger densities, the pairing acquires 
contributions from multiple Fermi points. Similarly, in 
higher dimensions and at low density, the problem has 
rotational symmetry, and the FFLO order parameter can 
be faithfully described by the modulus of a single wave 
vector |Q|. However, at higher densities the Fermi sur- 
face(s) acquires a diamond-like shape, and a multi-modal 
description is expected to be more accurate [26]. This 
would translate into complex real-space modulations of 
the order parameter, such as those observed in Ref. [13]. 
Beyond ID, nesting between bands with opposite spin is 
weak [5] , particularly at high densities and high magnetic 
fields, when the Fermi surfaces have very different shapes, 
with a large mismatch in Fermi velocities. Thus, the 
FFLO state becomes unfavorable compared to a normal 
polarized state. The ladder system, however, is strongly 
nested to the effect that we find a pairing instability of 
the FFLO-type at all p and n studied, with a much richer 
structure than in the single-chain case. 
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